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Abstract. We show that the effect of principal pivot transform on the 
nullity values of the principal submatrices of a given (square) matrix is 
described by the symmetric difference operator (for sets). We consider 
its consequences for graphs, and in particular generalize the recursive 
relation of the interlace polynomial and simplify its proof. 



1 Introduction 

Principal pivot transform (PPT, or simply pivot) is a matrix transformation 
operation capable of partially (component-wise) inverting a given matrix. PPT 
is originally motivated by the well-known linear complementarity problem [20], 
and is applied in many other settings such as mathematical programming and 
numerical analysis, see [19] for an overview. 

A natural restriction of pivot is to graphs (with possibly loops), i.e., sym- 
metric matrices over F2. For graphs, each pivot operation can be decomposed 
into a sequence of elementary pivots. There are two types of elementary pivot 
operations, (frequently) called local complementation and edge complementa- 
tion. These two graph operations are also (in fact, originally) defined for simple 
graphs. The operations are similar for graphs and simple graphs, however, for 
simple graphs, applicability is less restrictive. Local and edge complementation 
for simple graphs, introduced in [16] and [5] respectively, were originally mo- 
tivated by the study of Euler circuits in 4-regular graphs and by the study of 
circle graphs (also called overlap graphs) as they model natural transformations 
of the underlying circle segments. Many other applications domains for these op- 
erations have since appeared, e.g., quantum computing [21], the formal theory 
of gene assembly in ciliates [11] (a research area within computational biology), 
and the study of interlace polynomials, initiated in [1]. In many contexts where 
local and edge complementation have been used, PPT has originally appeared 
in disguise (we briefly discuss some examples in the paper). 

In this paper we show that the pivot operator on matrices A (over some 
field) and the symmetric difference operator on sets Y have an equivalent effect 
w.r.t. the nullity value of the principal submatrices A[Y] of A. We subsequently 
show that this nullity invariant can be formulated in terms of (a sequence of) 
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set systems. Furthermore we discuss its consequences for pivot on graphs and 
in particular apply it to the interlace polynomial. It was shown in [3] that the 
interlace polynomial, which is defined for graphs, fulfills a characteristic recur- 
sive relation. We generalize the notion of interlace polynomial and its recursive 
relation to square matrices in general. In this way, we simplify the proof of the 
(original) recursive relation for interlace polynomials of graphs. Also, in Sec- 
tion 3, we recall a motivation of pivot applied to overlap graphs, and relate it to 
the nullity invariant. 

2 Notation and Terminology 

A set system (over V) is a tuple M = (V, D) with V a finite set, called the 
domain of M, and D C 2 V a family of subsets of V. To simplify notation we 
often write X £ M to denote X £ D. Moreover, we often simply write V to 
denote the domain of the set system under consideration. We denote by © the 
logical exclusive-or (i.e., addition in F 2 ), and we carry this operator over to sets: 
for sets A,B CV, A® B is the set defined by x £ A © B iff (x £ A) © (x £ B) 
for x £ V. For sets, the © operator is called symmetric difference. 

We consider matrices and vectors indexed by a finite set V. For a vector v 
indexed by V, we denote the element of v corresponding to i £ V by v[i}. Also, 
we denote the nullity (dimension of the null space) and the determinant of a 
matrix A by n(A) and det(A) respectively. For X CV, the principal submatrix 
of A w.r.t. X is denoted by A[X\. 

We consider undirected graphs without parallel edges, however we do allow 
loops. Hence a graph G = (V, E) can be considered a symmetric V x F-matrix 
A = (a UiV ) over F 2 (the field having two elements): for u £ V, {u} £ E (i.e., u 
has a loop in G) iff a u , u — 1, and for u, v £ V with a/t, {u, v} £ E iff a UjV = 1. 
We denote the set of edges of G by E(G). We often make no distinction between 
G and its matrix representation A. Thus, e.g., we write n{G) = n(A), and, for 
X C V, G[X] = A [X], which consequently is the subgraph of G induced by X. 
Note that as G is represented by a matrix A over F 2 , n(G) is computed over F 2 . 
Also, for Y C V, we define G\Y = G[V \Y]. In case Y = {v} is a singleton, 
to simplify notation, we also write G\Y = G\v. Similar as for set systems, we 
often write V to denote the vertex set of the graph under consideration. 

3 Background: Nullity and Counting Closed Walks 

In this section we briefly and informally discuss an application of principal pivot 
transform where nullity plays an important role. In [9] a first connection be- 
tween counting cycles and the nullity of a suitable matrix was established. It is 
shown in that paper that the number of cycles obtained as the result of applying 
disjoint transpositions to a cyclic permutation is described by the nullity of a 
corresponding "interlace matrix" . 

It has been recognized in [18] that the result of [9] has an interpretation 
in terms of 2-in, 2-out digraphs (i.e., directed graphs with 2 incoming and 2 
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Fig. 1. The overlap graph of s = 146543625123. 




Fig. 2. A 2-in, 2-out digraph. 



outgoing edges for each vertex), linking it to the interlace polynomial [2]. We 
discuss now this interpretation in terms of 2-in, 2-out digraphs and subsequently 
show the connection to the pivot operation. 

Let V = {1, 2, 3, 4, 5, 6} be an alphabet and let s = 146543625123 be a double 
occurrence string (i.e., each letter of the string occurs precisely twice) over V. 
The overlap graph O s corresponding to s has V as the set of vertices and an 
edge {u, v} precisely when u and v overlap: the vertices u and v appear cither 
in order u,v,u,v or in order v,u,v,u in s. The overlap graph O s is given in 
Figure 1. One may verify that the nullity of O s is n(O s ) — 0. Consider now the 
subgraph O s [X] of O s induced by X = {3,4,5,6}. Then it can be verified that 
n(O s [X})=2. 

We discuss now the link with 2-in, 2-out digraphs (only in this section we 
consider digraphs). Let G be the 2-in, 2-out digraph of Figure 2 with V = 
{1,2, 3, 4, 5, 6} as the set of vertices. Although our example graph does not have 
parallel edges, there is no objection to consider such "2-in, 2-out multidigraphs" . 
Notice that the double occurrence string s — 146543625123 considered earlier 
corresponds to an Eulcr circuit C of G. We now consider partitions P of the edges 
of G into closed walks (i.e., cycles where repeated vertices are allowed). Note that 
there are 2^1 such partitions: if in a walk passing through vertex v we go from 
incoming edge e of v to outgoing edge e' of v, then necessarily we also walk in 
P from the other incoming edge of v to the other outgoing edge of v. Hence 
for each vertex there are two "routes". Let P now be the the partition of the 
edges of G into 3 closed walks as indicated by Figure 3 using three types of line 
thicknesses. Then P follows the same route as the Euler circuit (corresponding 
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Fig. 3. Partition of the edges of a 2-in, 2-out digraph into three closed walks. 



to) s in vertices {1, 2}, while in the other vertices X = {3, 4, 5, 6} it follows the 
other route. We say that P is induced by X in s. 

Theorem 1 in [9] now states (applying it to the context of 2-in, 2-out digraphs) 
that the number of closed walks of a partition P of edges induced by X in s is 
n(O s [X]) + 1. In our case we have indeed \P\ = 3 and n(O s [X]) = 2. 

The pivot operation, which is recalled in the next section, has the property 
that it can map Sl into S2 for any two double occurrence strings s\ and 
S2 that correspond to Euler circuits of a 2-in, 2-out digraph G, see, e.g., the 
survey section of [6]. For example, the partition of edges induced by {1,3} in 
s corresponds to a single closed walk which may be described by the double 
occurrence string s' = 123625146543. It then holds that O s i is obtained from O s 
by pivot on {1,3}, denoted by O s < — O s * {1,3}. We notice that the partition 
induced by {1,3} ® {3,4,5,6} = {1,4,5,6} in s' is equal to the partition P 
induced by {3,4,5,6} in s depicted in Figure 3. Hence n(O s * Y[Y © X}) = 
n(O s [X]) for X = {3, 4, 5, 6} and Y = {1, 3}. In Theorem 5 below we prove this 
property for arbitrary X and Y and for arbitrary square matrices (over some 
field) instead of restricting to overlap graphs O s . 

4 Pivot 

In this section we recall principal pivot transform (pivot for short) for square 
matrices over an arbitrary field in general, see also [19]. 

Let A be a V x V- matrix (over an arbitrary field) , and let X C V be such that 
the corresponding principal submatrix A[X] is nonsingular, i.e., dct A[X] 0. 
The pivot of A on X, denoted by A * X, is defined as follows. If P = A[X] and 
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Q is called the Schur complement of P in A. 
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The pivot can be considered a partial inverse, as A and A*X are related by 
the following equality, where the vectors x\ and y\ correspond to the elements of 
X. This equality is characteristic as it is sufficient to define the pivot operation, 
see [19, Theorem 3.1]. 

Note that if det A j£ 0, then A * V = A^ 1 . Also note by Equation (2) that the 
pivot operation is an involution (operation of order 2), and more generally, if 
(A * X) * Y is defined, then it is equal to A * (X © Y). 



5 Nullity Invariant 

It is well known that any Schur complement in a matrix A has the same nullity as 
A itself, see, e.g., [22, Section 6.0.1]. See moreover [22, Chapter 0] for a detailed 
historical account of the Schur complement. We can rephrase the nullity property 
of the Schur complement in terms of pivot as follows. 

Proposition 1 (Nullity of Schur complement). Let A be a V x V -matrix, 
and let X CV such that A[X] is nonsingular. Then n(A* X[V\X)) = n(A[V}). 

The following result is known from [20] (see also [10, Theorem 4.1.2]). 

Proposition 2. Let A be a V x V -matrix, and let X CV be such that A[X] is 
nonsingular. Then, for Y CV, 

det(A*X)[Y] =detA[X ®Y]/detA[X]. 

As a consequence of Proposition 2 we have the following result. 

Corollary 3. Let A be a V x V -matrix, and let X C V be such that A[X] 
is nonsingular. Then, for Y C V, (A * X)[Y) is nonsingular iff A[X © Y] is 
nonsingular. 

We will now combine and generalize Proposition 1 and Corollary 3 to obtain 
Theorem 5 below. 

We denote by A\X the matrix obtained from A by replacing every row t>J of 
A belonging to x G V \ X by i^ where i x is the vector having value 1 at element 
x and elsewhere. 

Lemma 4. Let A be a V x V -matrix and X CV. Then n{A$X) = n(A[X}). 

Q 



Proof. By rearranging the elements of V, A is of the following form 

Q 



where A[X] = P. Now A$X is 



R 

where / is the identity matrix of suitable 



size. We have n(P) = n(A$X). □ 
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We are now ready to prove the following result, which we refer to as the 
nullity invariant. 



Theorem 5. Let A be a V x V -matrix, and let X C V be such that A[X] is 
nonsingular. Then, for Y C V, n((A * X){Y]) = n(A[X © Y}). 

Proof. We follow the same line of reasoning as the proof of Parsons [17] of Propo- 
sition 2 (see also [10, Theorem 4.1.1]). Let Ax = y. Then 



((AiX)xM = ■ 

As, by Equation (2), 

((A * X)(A$X)x)[i] 



y[i] itieX, 
x\i] otherwise. 



x[i] if i S X, 
y[i] otherwise, 



we have, by considering each of the four cases depending on whether or not i in 
X and i in Y separately, 



(((A*XW)(A$X)x)\i] 



y[i] iiieX®Y, 
x[i] otherwise. 



Thus we have ((A * X)$Y)(ApC) = A$(X © Y). By Lemma 4, n{A$X) = 
n{A[X\) = 0, and therefore A\X is invertible. Therefore n((A*X)$Y) = n(A$(X® 
Y)), and the result follows by Lemma 4. □ 

By Theorem 5, we see that the pivot operator *I on matrices and the sym- 
metric difference operator ®X on sets have an equivalent effect on the nullity 
values of principal submatrices. 

Note that Theorem 5 generalizes Corollary 3 as a matrix is nonsingular iff 
the nullity of that matrix is (the empty matrix is nonsingular by convention) . 
One can immediately see that Theorem 5 generalizes Proposition 1. 

Also note that by replacing Y by V \ Y in Theorem 5, we also have, equiva- 
lent^, n((A * X)[X © Y}) = n(A[Y]). 

The "Nullity Theorem" [13, Theorem 2], restricted to square principal sub- 
matrices, states that if A is an invertible V x T^-matrix, then, for Y C V, 
n(^ _1 [F]) = n(A[V\Y]). Note that this is implied by Theorem 5 as A*V = A' 1 . 

/125\ 

Example 6. Let V — {a, b, c} and let A be the V x ^-matrix 14 2 over 

\321/ 

Q where the columns and rows are indexed by a, 6, c respectively. We see that 
/4 2\ 

A[{b, c}] — ( 2 x ) anc ^ therefore n(A[{b, c}]) = 1. Moreover, for X = {a, b}, the 
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Fig. 4. Graphs G and G * X of Example 7. 



columns of A[X] are independent and thus det A[X] ^ 0. We have therefore that 
A * X is defined, and it is given below. 



A*X 



By Theorem 5, we have n(A[{b, c}]) = * X[X © {&, c}]) = n(A * X[{a, c}]). 
Therefore n(^4 * JT[{a, c}]) = 1. This can easily be verified given A * X[{a, c}] = 

2 -8 

5 -20 




□ 



It is easy to verify from the definition of pivot that A * X is skew-symmetric 
whenever A is. In particular, if G is a graph (i.e., a symmetric matrix over F 2 ), 
then G * X is also a graph. For graphs, all matrix computations, including the 
determinant, will be over F 2 . 

Example 7. Let G be the graph given on the left-hand side of Figure 4. Let 

/001\ 

X = {I, 2, 3}. Then the X x X-matrix belonging to G[X] is 1 1 where the 

Vlll/ 

columns and rows represent vertices 1, 2, 3, respectively. We see that the columns 
of G[X] are independent (over F 2 ) and therefore detG[X] = 1. Consequently 
G * X is defined and the graph is given on the right-hand side of Figure 4. 
Let now Y = {1,4}. We see that G[Y] is a discrete graph (i.e., the graph has 
no edges). Therefore n(G[Y}) = 2. Now by Theorem 5, we have n(G[Y}) = 
n(G * X[X © Y}) = n(G * X{{2, 3, 4}]). One may verify that removing vertex 1 
from G * X indeed obtains a graph of nullity 2. □ 



6 Set Systems 

Let A be a V x ^-matrix. Let Ma = {V,D) be the set system with X e D 
iff A[X] is nonsingular. Set system A4a turns out to fulfill a specific exchange 
axiom if A is (skew-)symmetric, making it in this case a delta-matroid [4] (we 
will not recall its definition here as we do not use this notion explicitly). 
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Let M = (V, D) be a set system. We define for XCV, the pivot (often called 
twist) ofMonl, denoted M*X, by (V, D*X) where D*X = {Y®X \ Y G M}. 
By Corollary 3 it is easy to verify, see [14], that the operations of pivot on set 
systems and matrices match, i.e., Ma * X = Ma*x if the right-hand side is 
defined (i.e., if X G Ma)- 

Theorem 5 allows now for a generalization of this result from the set system 
Ma oi nullity to a "sequence of set systems" Va for each possible nullity i. 
We formalize this as follows. 

For a finite set V, we call a sequence P = (Pq,P\, . . . , P n ) with n — \V\ and 
Pi C V for alH G {0, . . . , n} a partition sequence (over V) if the nonempty Pi's 
form a partition of 2 V . Regarding P as a vector indexed by {0, ... , n}, we denote 
Pi by P[i]. Moreover, we define for partition sequence P and X C V, the pw>o£ of 
P on X, denoted by P*X, to be the partition sequence (Pq*X, Pi*X 7 . . . , P n *X). 
Also, we call the vector (|P |, |Pi|, . . . , |P„|) of dimension n + 1, denoted by ||P||, 
the norm of P. Clearly, ||P|| = ||P * X\\, i.e., the norm of P is invariant under 
pivot. 

For a Y x Y-matrix A we denote by Va the partition sequence over V where 
X e Pa[i\ iff n(A [X]) = i. As nullity corresponds to a non-zero determinant 
(this holds also for as det A[0] = 1 by convention), we have Ma = {V 1 P^O]). 

We now have the following consequence of Theorem 5. Note that X 6 Va[0] 
iff A * X is defined. 

Theorem 8. Let A be a V x V -matrix, and X G V A [0]- Then V A *x =Va*X. 

Proof. By Theorem 5 we have for all i € {0, . . . , n}, Y € Pa*x[«] iff * 
X)[F]) = i iff n(A[A © Y]) = i iff X © Y e P A \i} iff Y e P A [i] * X. □ 

Since the norm of a partition sequence is invariant under pivot, we have by 
Theorem 8, |[Pa|| = ||Pa*x||- Therefore, for each i G {0, ...,n}, the number 
of principal submatrices of A of nullity i is equal to the number of principal 
submatrices of A * X of nullity i. 

For X C V, it it is easy to see that Va[x] 1S obtained from Va by removing 
all Y G Va[i] containing at least one element outside X: VA[x]{i] = {Z C X \ 
Z €V A [i\}{or &\ii€{0,...,\X\}. 

Example 9. For matrix A from Example 6, we have Va = (Po, Pi, P2, Ps) with 
P = 2 y \ {{fe, c}}, Px = {{6, c}}, and P 2 = P 3 = 0. □ 

Example 10. For graph G from Example 7, depicted on the left-hand side of 
Figure 4, we have V G = (P , Pi, P 2 , P 3 , P 4 ) with 

P O = {0,{2},{3},{1,3},{2,4},{3,4},{1,2,3},{1,2,3,4}}, 
Pi = {{1},{4},{1,2},{2,3},{1,2,4},{1,3,4},{2,3,4}}, 

P 2 = {{1,4}},P 3 = P4 = 0. 
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By Theorem 8 we have for G * X with X = {1, 2, 3}, depicted on the right-hand 
side of Figure 4, V G *x = {P^P^P^, Pi) where 

i^ = {0, {2}, {4}, {1,2}, {1,3}, {1,2,3}, {1,2,4}, {1,3,4}}, 
P{ = {{1}, {3}, {1, 4}, {2, 3}, {2, 4}, {3, 4}, {1, 2, 3, 4}}, 
P^ = {{2,3,4}}, P^ = Pi = 0. 

We have p> G || - \\V G *x\\ = (8,7,1,0,0). □ 

Example 11. In the context of Section 3, where matrix A an overlap graph O s 
for some double occurrence string s, we have that HPoJIH is the number of 
partitions of the edges of the 2-in, 2-out digraph D corresponding to s into 
closed walks of D, such that the number of closed walks is precisely i + 1. The 
value HPoJIlO] is therefore the number of Euler circuits in D. □ 

7 Elementary Pivots on Graphs 

From now on we consider pivot on graphs (i.e., symmetric V x ^-matrices over 
F2), and thus on all matrix computations will be over F2. Hence for graph G, 
Mg = (V,D G ) is the set system with X £ D G iff det(G[X]) = 1. Also, G can 
be (re) constructed given M G . Indeed, {u} is a loop in G iff {u} £ D G , and 
{u,v} is an edge in G iff {{u,v} £ D G ) © (({«} 6 D G ) A ({v} £ D G )), see [7, 
Property 3.1]. Therefore, the function M.(.) assigning to each graph G the set 
system M G is an injective function from the family of graphs to the family of 
set systems. It this way the family of graphs may be regarded as a subclass of 
the family of set systems. Note that Ai G * X is defined for all X C V, while 
pivot on graphs G*X is defined only if X £ M G (or equivalently, G M G *X). 

In this section we recall from [14] that the pivot operation on graphs can 
be defined as compositions of two graph operations: local complementation and 
edge complementation. 

The pivots G*X where X is a minimal element of M G \{0} w.r.t. inclusion 
are called elementary. It is noted in [14] that an elementary pivot X on graphs 
corresponds to either a loop, X = {u} e E(G), or to an edge, X — {u,v} e 
E(G), where both vertices u and v are non-loops. Thus for Y G M G , if G[Y] has 
elementary pivot Xi, then Y\X\ = Y®X\ £ M G *x x ■ In this way, each Y £ M. G 
can be partitioned Y = X x U • • • U X n such that G *Y = G * (Xi • • • ® X n ) = 
(• • • (G * X\) ■ ■ ■ * X n ) is a composition of elementary pivots. Consequently, a 
direct definition of the elementary pivots on graphs G is sufficient to define the 
(general) pivot operation on graphs. 

The elementary pivot G*{u} on a loop {u} is called local complementation. It 
is the graph obtained from G by complementing the edges in the neighbourhood 
N G (u) — {v £ V I {u,v} £ E(G),u ^ v} of u in G: for each v,w £ N G (u), 
{v, w} £ E(G) iff {v, w} <£ E(G * {«}), and {v} £ E{G) iff {v} E(G * {u}) 
(the case v = w). The other edges are left unchanged. 

The elementary pivot G * {it, v} on an edge {u, v} between distinct non-loop 
vertices u and v is called edge complementation. For a vertex x consider its closed 
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Fig. 5. Pivoting on an edge {u, v} in a graph with both u and v non loops. Connection 
{x, y} is toggled iff x € Vi and y G Vj with i / j. Note u and t> are connected to all 
vertices in Vz, these edges are omitted in the diagram. The operation does not affect 
edges adjacent to vertices outside the sets Vi,V2,V3, nor does it change any of the 
loops. 




Fig. 6. The orbit of a graph under pivot. Only the elementary pivots are shown. 

neighbourhood Nq(x) = Ng(x) U {x}. The edge {u, v} partitions the vertices of 
G connected to u or v into three sets V\ — N' G (u) \ Nq(v) 7 V2 — Nq(v) \ Nq(u), 
V 3 = N^(u) n N^(v). Note that u, v e V 3 . 

The graph G * {u, v} is constructed by "toggling" all edges between different 
Vi and Vj-. for {x,y} with x 6 Vi, y G Vj and i ^ j: {x,y} G E{G) iff {x,y} ^ 
E(G[{u, v}]), see Figure 5. The other edges remain unchanged. Note that, as a 
result of this operation, the neighbours of u and v are interchanged. 

Example 12. Figure 6 depicts an orbit of graphs under pivot. The figure also 
shows the applicable elementary pivots (i.e., local and edge complementation) 
of the graphs within the orbit. □ 

fnterestingly, in many contexts, principal pivot transform originally appeared 
in disguise. For example, PPT was recognized in [15] as the operation underlying 
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the recursive definition of the interlace polynomial, introduced in [1]. We will 
consider the interlace polynomial in the next section. Also, e.g., the graph model 
defined in [12] within the formal theory of (intramolecular) gene assembly in cil- 
iates turns out to be exactly the elementary pivots, as noted in [8]. Furthermore, 
the proof of the result from [9], connecting nullity to the number of cycles in 
permutations, as mentioned in Section 3, implicitly uses the Schur complement 
(which is an essential part of PPT) . 

8 The Interlace Polynomial 

The interlace polynomial is a graph polynomial introduced in [1,2]. We follow 
the terminology of [3]. The single- variable interlace polynomial (simply called 
interlace polynomial in [2]) for a graph G (with possibly loops) is defined by 

q (G) = J2(y-V n(G[s]) - 

scv 

It is is shown in [3] that the interlace polynomial fulfills an interesting recur- 
sive relation, cf. Proposition 15 below, which involves local and edge comple- 
mentation. As we consider here its generalization, principal pivot transform, it 
makes sense now to define the interlace polynomial for V x V-matrices (over 
some arbitrary field) in general. Therefore, we define the interlace polynomial 
for V x ^-matrix A as 

q (A) = J2(y-i) n{Als]) - 

scv 

We may (without loss of information) change variables y := y — 1 in the definition 
of the interlace polynomial to obtain 

q '(A) = J2y n(A[s]) - 

scv 

As q(A) (and q'(A)) deals with nullity values for (square) matrices in general, 
one can argue that the nullity polynomial is a more appropriate name for these 
polynomials. 

We see that the coefficient ai of term a,y* of q'(A) is equal to H^aHI*] (the 
element of \\Va\\ corresponding to i) for all i € {0, . . . , n}. Therefore, we have 
for matrices A and A', q(A) = q(A') iff q'(A) = q'(A') iff \\V A \\ = \\Va>\\- 

Example 13. Let O s be the overlap graph for some double occurrence string s, 
and let a* be the coefficient a* of term a^y* of q'(O s ). We have, see Example 11, 
that a, is equal to the number of partitions of the edges of the 2-in, 2-out digraph 
D corresponding to s into closed walks of D, such that the number of closed 
walks is precisely i + 1. More specifically, ao is the number of Euler circuits in 
D. The interlace polynomial is originally motivated by the computation of these 
coefficients at of 2-in, 2-out digraphs, see [2]. □ 
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It is shown in [2] that the interlace polynomial is invariant under edge com- 
plementation. By Theorem 8 we see directly that this holds for pivot in general: 
\\Pa*x\\ = \\V A \\ and cquivalcntly q(A * X) = q(A). 

Furthermore, we show that q(A) fulfills the following recursive relation. 

Theorem 14. Let A be aV x V -matrix (over some field), let X QV with A[X] 
nonsingular, and let u e X . We have q(A) = q(A \ u) + q(A * X \ u). 

Proof. Let Va = (Po> Pi, ■ ■ ■ , Pn)- Since X is nonempty and A[X] is nonsingular, 
P„ = 0. Let R = (P , Pi, ■ ■ ■ , Pn-i)- Let Z e Pi for i e {0, 1, . . . , n - 1}. We 
have Z <Z V does not appear in Va\ u iff u E Z iff u ^ Z Q) X iff Z Q) X does 
appear in Va*x\u- Hence ||P|| = ||7- > J 4\ tt || + ||Pa*x\m|| (point-wise addition of the 
two vectors), and the statement holds. □ 

The recursive relation for the single-variable interlace polynomial in [3] is 
now easily obtained from Theorem 14 by restricting to the case of elementary 
pivots on graphs. 1 

Proposition 15 ([3]). Let G be a graph. Then q(G) fulfills the following con- 
ditions. 

1. q(G) = q(G \ u) + q(G * {u, v} \ u) for edge {u, v} in G where both u and v 
do not have a loop, 

2. q{G) = q(G \ u) + q{G * {u} \u) if u has a loop in G, and 

3. q(G) =y n if G is a discrete graph with n vertices. 

Proof. Conditions (1) and (2) follow from Theorem 14 where A is a graph G, 
and X is an elementary pivot (i.e., X — {u} is a loop in G or X = {u,v} is an 
edge in G where both u and v do not have a loop, see Section 7). Finally, if G 
is a discrete graph with n vertices, then, for all Y C V, Y e P\y\- Consequently, 
|P| = (™). Thus, q'{G) = (y + 1)" and therefore q(G) =y n . □ 

9 Discussion 

We have shown that the pivot operator *X on matrices A and the symmetric 
difference operator (BX on sets Y have an equivalent effect w.r.t. the nullity value 
of the principal submatrices A[Y] of A. This nullity invariant may be described 
in terms of partition sequences Va, where the sets Y CV are arranged according 
to the nullity value of A[Y]. We notice that interlace polynomial of a graph G 
corresponds to the norm \\Vg\\ of the partition sequence of G (where G is con- 
sidered as a matrix). Hence we (may) naturally consider interlace polynomials 
for square matrices in general, and obtain a recursive relation for these gener- 
alized interlace polynomials. In this way, we simplify the proof of the (original) 
recursive relation for interlace polynomials of graphs. 

1 We use here the fact observed in [15] that the operations in the recursive relations of 
[3] are exactly the elementary pivots of Section 7, assuming that the neighbours of 
u and v are interchanged after applying the "pivot" operation of [3] on edge {u,v}. 
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